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ABSTRACT 

An N-body hybrid simulation, integrating both massive and tracer particles, of a 
Galactic disk is used to study the stellar phase space distribution or velocity distri- 
butions in different local neighborhoods. Pattern speeds identified in Fourier spec- 
trograms suggest that two-armed and three-armed spiral density waves, a bar and a 
lopsided motion are coupled in this simulation, with resonances of one pattern lying 
near resonances of other patterns. We construct radial and tangential (uv) velocity 
distributions from particles in different local neighborhoods. More than one clump 
is common in these local velocity distributions regardless of the position in the disk. 
Features in the velocity distribution observed at one galactic radius are also seen in 
nearby neighborhoods (at larger and smaller radii) but with shifted mean v values. 
This is expected if the v velocity component of a clump sets the mean orbital galac- 
tic radius of its stars. We find that gaps in the velocity distribution are associated 
with the radii of kinks or discontinuities in the spiral arms. These gaps also seem 
to be associated with Lindblad resonances with spiral density waves and so denote 
boundaries between different dominant patterns in the disk. We discuss implications 
for interpretations of the Milky Way disk based on local velocity distributions. Veloc- 
ity distributions created from regions just outside the bar's Outer Lindblad resonance 
and with the bar oriented at 45° from the Sun-Galactic center line more closely re- 
semble that seen in the solar neighborhood (triangular in shape at lower uv and with 
a Hercules like stream) when there is a strong nearby spiral arm, consistent with the 
observed Centaurus Arm tangent, just interior to the solar neighborhood. 



1 INTRODUCTION 



1999 



2006 



The velocity distribution of stars in the solar neighbor- 
hood contains structure that has been particularly clearly 



Fux|2001||Quillen fc Minchev|2005||Minchev fe Quillen 



Antoja et al.| |2009 Gar dner et al.|[2010[ |Lepine et al 



revealed from Hipparcos observations (Dehnen & Binney 
1998 Famaey et al.|2004 l. Much of this structure was previ- 



20101. Mergers and orbiting satellite subhalos can leave be- 
hind stellar streams (e.g., Bekki & Freeman 2003, Meza et 



al. 200m IHelmi et al. 2006 Gomez et al. 20101. In a local 



ously associated with moving groups ( Eggen| 19*96 1 . Moving 
groups are associations of stars that are kinematically sim- 
ilar or have similar space motions. Young, early-type stars 
can be moving together because they carry the kinematic 
signature of their birth ( Eggen 1996 1. However, a number of 
the kinematic clumps identified through the kinematic stud- 



ies also contain later type and older stars ( Dehnen & Binney 
T9981 I Famaey et al.|[2U58| |Nordstrom et al.||2004| |Arifyanto 



fc Fuchs|2006 l. 



In a small region or neighborhood of a galaxy, structure 
in the stellar velocity distribution can be caused by differ- 
ent dynamical processes. Bar or spiral patterns can induce 
structure in the phase space distribution, particularly in the 
vicinity of Lindblad resonances ( |Yuan fc Kuo|1997||Dehnen| 



neighborhood, these could be seen as increases in the stellar 
number density at a particular velocity. 

The velocity vector of a star in the solar neighborhood 
can be described in Galactocentric coordinates with compo- 
nents (it, v, w). Here u is the radial component with positive 
u for motion toward the Galactic center. The v component 
is the tangential velocity component subtracted by the ve- 
locity of a particle in a circular orbit computed from the 
azimuthally averaged midplane gravitational potential. The 
v component is positive for a particle with tangential com- 
ponent larger than the circular velocity that is moving in 
the same direction as disk rotation. The w component is the 
velocity in the z direction. 

IQuillen fc Minchev| ( |2005[ ) showed that Lindblad reso- 
nances associated with a spiral pattern, given by their ap- 
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proximate v velocity component, were likely locations for 
deficits in the velocity distribution in the solar neighborhood 
(also see Fux|2001 Lepine et al.|2010| ). On either side of the 



resonance the shape of orbits changes. For example closed 
or periodic orbits can shift orientation so they are parallel 
to a bar on one side of the resonance and perpendicular on 
the other side. On resonance there may be no nearly circu- 
lar orbits. As a consequence in a specific local neighborhood 
there can be a particular velocity associated with resonance 
that corresponds only to orbits with high epicyclic varia- 
tions. These may not be well populated by stars leading to 
a deficit of stars in the local velocity distribution at this 
particular velocity. This picture is also consistent with the 
division between the Hercules stream and the thin disk stars 
that is likely caused by the m = 2 outer Lindblad resonance 
(OLR) with the Galactic b ar ( |Dehnen|2000] |Minchev et al 



2007| |Gardner et al.||2010 l 



Structure seen in the velocity distribution can also be 
caused by recent large scale perturbations. Perturbations 
on the Galactic disk caused by bar formation or a galactic 
encounter, can perturb the action angle distribution. Subse- 
quent phase wrapping can induce structure in the velocity 
distribution ( |Minchev et al.|2009|[Quillen et al.|2009[ ). Tran- 



sient spiral density waves can also leave behind structure in 
the local phase space distribution ( De Simone et aL||2004 | 



The local velocity distributions can be particularly complex 
when both spiral and bar perturbations are present or when 



multiple spiral patterns are present l 


Quillen 2003 
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illen||2006 |Chakrabarty & Sideris 2008 


Antoja et al. 


Minchev & Famaey|2010 
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2011). 



The variety of explanations for structure in the regional 
phase space distribution presents a challenge for interpreta- 
tion of current and forthcoming radial velocity and proper 
motion surveys. We may better understand the imprint of 
different dynamical structures in the phase space distribu- 
tion by probing simplified numerical models. Most numerical 
studies of structure in localized regions of the galactic disk 
phase space distribution has been done using test particle in- 
tegrations (e.g., Dehnen 2000 De Simone et al. 2004 Quillen 



2009 



fc Minchev|2005| |Chakrabarty fc Sideris|2008| |Antoja et al 



Minchev et aL|2009 Quillen et aL|2009 Gardner et al 



20101. Few studies have probed the local phase space distri- 
bution in an N-body simulation of a galactic disk (but see 
Fux|2000 1 because of the large number of particles required 
to resolve structure both in physical and velocity space. In 
this paper we use an N-body hybrid simulation, integrat- 
ing both massive and tracer particles, that exhibits both 
spiral and bar structures to probe the velocity distribution 
in different regions of a simulated disk galaxy. Our goal is 
to relate structure in the local phase space distribution to 
bar and spiral patterns that we can directly measure in the 
simulation. 



2 SIMULATIONS 

In this section we describe our N-body hybrid simulations 
and their initial conditions. We also describe our procedure 
for making spectrograms of the Fourier components so that 
we can identify pattern speeds in the simulations. 



2.1 The N-body hybrid code 

The N-body integrator used is a direct-summation code 
called 0-GRAPE dHarfst et al. 1120071) that employs a 4-th 



order Hermite integration scheme with hierarchical com- 
mensurate block times steps ( |Makino fc Aa rscth 1992JI. In- 
stead of using special purpose GRAPE hardware we use 
the Sapporo subroutine library ( jGaburov et al.||2009[) that 
closely matches the GRAPE-6 subroutine library (Makino 
et al.|2 003 ) but allows the force computations to be done on 
graphics processing units (GPUs). The integrator has been 
modified so that massless tracer particles can simultaneously 



be integrated along with the massive particles (Comparetta 
|fc Qui llen 2010|. The tracer particles help us resolve struc- 
ture in phase space without compromising the capability to 
carry out a self-consistent N-body simulation. The simula- 
tions were run on a node with 2 NVIDIA Quadro FX 5800 
GPUs. The cards are capable of performing double precision 
floating point computations, however the computations were 
optimized on these GPUs to achieve double precision accu- 
racy with single precision computations using a corrector 
( |Gaburov et al.|2009[ ). The Quadro FX cards were designed 
to target the Computer-Aided Design and Digital Content 
Creation audience and so were slower but more accurate 
than other similar video cards even though they did not 
impliment error correction code. 



2.2 Initial Conditions for the Integration 

The initial conditions for a model Milky Way galaxy were 
made with numerical phase phase distribution functions us- 



ing the method discussed by Widrow et al. ( 2008 1 and their 
numerical routines which are describ ed by |Widrow et al. 
|2008| ; |Kuijken fc Dubinski| ( |1995[ ); |Widrow fc Dubinski 
(12005 1. This code, called GalacticICS, computes a gravita- 



tional potential for bulge, disk and halo components, then 
computes the distribution function for each component. Par- 
ticle initial conditions are then computed for each compo- 
nent. The galactic bulge is consistent with a Sercic law for 
the projected density. The halo density profile is described 
by 5 parameters, a halo scale length, a cusp exponent, a ve- 
locity dispersion, and an inner truncation radius and trunca- 
tion smoothness, so it is more general than an NFW profile. 
The disk falls off exponentially with radius and as a sech 2 
with vertical height. We have adopted the parameters for the 



Milky Way model listed in Table 2 by [Widrow et al.| ( |2008 ). 
A suite of observational constraints are used to constrain 
the parameters of this model, including the Oort constants, 
the bulge dispersion, the local velocity ellipsoid, the outer 
rotation curve, vertical forces above the disk and the disk 
surface brightness profile. 

The number of massive particles we simulated for the 
halo, bulge and disk are 50,000, 150,000 and 800,000, re- 
spectively. The halo is live. The number of test (massless) 
particles in the disk is 3 million and exceeds the number of 
massive particles by a factor of 3. The total mass in disk, 
bulge and halo is 5.3 x 10 10 , 8.3 x 10 9 , and 4.6 x 10 11 M Q , 
respectively. The mass of halo, bulge and disk particles is 

9.2 x 10 6 ,5.5 x 10 4 , and 6.6 x 1O 4 M , respectively. Snap- 
shots were output every 5 Myrs and our simulation ran for 

1.3 Gyr. 

The smoothing length is 10 pc and is similar to the mean 
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disk interparticle spacing. We have checked that the center 
of mass of the simulation drifts a distance shorter than the 
smoothing length during the entire simulation. The parame- 
ters used to choose timesteps in the Hermite integrator were 
i) = 0.1 and r) a = 0.01 (see |Makino fc Aarseth||1992| [Harfat 
etl[L][2007| . The minimum timestep is proportional to y/rj 
and is estimated from the ratio of the acceleration to the 
jerk. Our value places the minimum possible timestep near 
the boundary of recommended practical The total energy 
drops by 0.2% during the simulation, at the boundary of 
recommended practice t hat suggests a maximum en ergy er- 



^^/E 1 < 1/y/N (see |Quinlan fc Tremaine|l992 



on how 

the shadow distance depends on noise and softening). 

In our goal to well resolve the galactic disk we have pur- 
posely undersampled both halo and bulge and so have in- 
troduced spurious numerical noise into the simulation from 
these populations. The product of the mass of halo particles 
and the halo density density exceeds that of the same prod- 
uct for the disk implying that numerical heating from the 
massive particles in the disk is lower than that induced by 
the live halo particles. The large timesteps, relatively small 
smoothing length and large halo particles imply that this 
simulation is noisier than many discussed in the literature. 

We first discuss the morphology of the simulated galaxy 
and then use the spectrograms constructed from the disk 
density as a function of time to identify patterns. 



3 MORPHOLOGY 

In Figures [I] we show the disk density projected onto the 
midplane in Cartesian coordinates for a few of the simula- 
tion outputs (snapshots). Disk particles are used to create 
these density histograms, but bulge and halo particles are 
neglected. In Figures[l]a,b both massive and tracer particles 
are shown. Figure [Tj; is equivalent to Figure [Tja except only 
massive disk particles are shown. We find no significant dif- 
ferences between the distribution of massive and tracer disk 
particles, consistent with the previous study using the same 
hybrid code ( Comparetta fc Quillen|2010| ). 



The numbers of massive particles is sufficiently low that 
numerical noise is present and amplified to become spiral 
structure. During the initial 40 snapshots (over a total time 
of 200 Myr) (shown in Figure [l^) two-armed and three- 
armed spiral armed waves grow in the region with radius 
5-12 kpc. At around a time T ~ 225 Myr an open and 
somewhat lopsided axisymmetric structure grows within a 
radius of 4 kpc that becomes the bar. At later times both 
two-armed and three-armed structure exists outside the bar 
(shown in Figure [Tja). 

As our halo is live (and not fixed) lopsided motions 
are not prevented in the simulation. Lopsided motions and 
three-armed structure has previously been observed in N- 
body simulations of isolated disks (e.g, see the simulations 
by Chilingarian et al. 2010 1. The bulge does not remain fixed 



during the simulation as a lopsided perturbation develops 
that moves its center slowly up to a distance of 0.5 kpc 



N-BODY Simulation Techniques: There's e 
right way and a wrong way, by Katz, N. 
http: / /supernova. lbl.gov/~evlinder/umass / sumold / com81x.ps 
" http://www.ifa.hawaii.edu/~barnes/ast626_09/NBody.pdf 



away from the initial origin. This distance is larger than 
the distance our center of mass drifts during the simulation 
(less than 10 pc). The lopsided motion may be related to 
the development of three-armed structures in the disk and 
asymmetries in the bar. Lopsidedness is common in spiral 



galaxies (see Jog & Combes 2009 for a review), and can 
spontaneously grow in N-body simulations of disks (Revaz 
fc Pfenniger|2004| [Saha et al.|2007[ ). If the pattern speed is 



slow, a lopsided mode can be long lasting ( |Ideta|2002[ ). There 
is no single strong frequency or amplitude associated with 
the lopsided motion in our simulation, though the motion 
is primarily counter-clockwise and moving with the direc- 
tion of galactic rotation. We have checked that the center of 
mass of the simulation drifts less than the smoothing length 
of 10 pc so the lopsided motion is not due to an unphysical 
drift in the total momentum in the simulation. However, the 
growth of the lopsided motion could be due to noise associ- 
ated with our small smoothing length and large halo mass 
particles. Simulations of a disk residing in a rigid halo would 
not exhibit lopsided motions. However as most galaxies dis- 
play lopsidedness we can consider it a realistic characteristic 
of a galactic disk. We will discuss this motion later on when 
we discuss the m = 1 Fourier component spectrograms. 

We also show the disk density in polar coordinates in 
Figures [2] and [3] These fi gures show differential density his- 
tograms as a function of radius and azimuthal angle (r, <j>) 
(or in cylindrical coordinates again projected onto the mid- 
plane). The bins to create these two plots are logarithmically 
spaced in radius. The disk density E(r, <j>) in these two plots 
has been normalized at each radius. These histogram show 
(E(r, ^)-E(r))__(r) _1 where E(r) is the mean density at ra- 
dius r averaged over the angle (f>. A logarithmic spiral would 
give high densities along a straight line in these plots, with 
slope set by the spiral arm pitch angle. Two armed struc- 
tures would correspond to two linear but sloped features, 
each separated by 180°. Trailing spiral structures have neg- 
ative slopes. A bar on these diagrams would correspond to 
two vertical features separated by 180° . The centroid of the 
galaxy bulge is subtracted prior to computing the density 
histograms shown in[2]and[3]in polar coordinates but not in 
Figure [l] that is in Cartesian coordinates. The times of the 
snapshots in Myr are shown in the upper right hand corners 
of each panel shown in Figures [I] - [3] 

In Figure [2] the differential surface density in polar co- 
ordinates is shown at three different times early in the simu- 
lation and during the early phase of spiral arm growth. The 
nearly vertical strip in the leftmost panel tilts and moves 
to the right at later times. In the right most panel it is 
a stronger perturbation and with a larger, trailing pitch 
angle. The spiral arm shown in the center of the leftmost 
panel becomes more tightly wound, and stronger (as seen 
on the right) as would be expected from a swing amplifi- 
cation mechanism. After about 50 Myr the spiral patterns 
are more slowly evolving. Other than during bar formation 
at about 330 Myr, and in the beginning of the simulation 
before 50 Myr, the spiral pitch angles do not seem to vary 
significantly. 

In Figure [3] the differential surface density is also shown 
at three different times but after bar formation. On the lower 
parts of the panels the bar is seen as pair of vertical strips. 
We can see that the bar length increases later in the simula- 
tion. Exterior to the bar, both two-armed and three-armed 
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patterns are seen and their patterns move more slowly than 
the bar. The winding or pitch angle of spiral features is about 
24° and a = d(/>/d\nr ~ 2.3. This pitch angle exceeds many 
but not all estimates for spiral arm pitch angles near the 
solar neighborhood (see Table 1 by Vallee |2008 most pitch 
angle estimates are ~ 13°). Other than during bar forma- 
tion and before 50 Myr, the spiral pitch angles do not seem 
to vary significantly. 

When animated the two-armed and three-armed fea- 
tures appear to constructively and destructively add with 
density peaks occurring at times and positions when two 
patterns lie on top of one another (e.g., as discussed by 
Henry et al.|2003||MeTdt et al.|2009| . The sum of two differ- 
ent patterns can cause individual arm peaks to shift position 
and temporarily appear more or less tightly wound. 

Kinks or changes in pitch angle of the spiral arms seen 
in the density images in Cartesian coordinates (in FigurefTJs) 
correspond to places where the spiral arms are discontinuous 
or change slope in the density distribution as seen in the 
polar plots (e.g., Figure |3j. For example, the spiral pattern 
just outside the bar has a lower pattern speed than the bar. 
Consequently the phase of the spiral arms varies with respect 
to the bar. At times spiral arms appear to be linked with 
the bar and at other times they appear separated from the 
bar. Spiral arms in the outer galaxy have slower pattern 
speeds than those in the inner galaxy. At times the spiral 
arms appear to be connected and at other times, armlets or 
kinks are seen. 



4 SPIRAL AND BAR PATTERNS IN THE 
SIMULATION 

Before we discuss velocity distributions extracted from local 
neighborhoods we first discuss patterns of waves measured 
from the simulations. 



4.1 Constructing Spectrograms 

We constructed spectrograms using the procedure outlined 



by Meidt et al. (20081 in their section 3.2 (also see section 
3.4 by Masset & Tagger 19971. In each radial bin, defined 



by a binning function b(r), and at each snapshot time t we 
compute 



Wc{r,t) = ^cos(m0<)6(r) 

i 

W?(r,t) = 5^sin(mfli)6(r) 



(1) 



where the sum is over all disk particles in the radial bin (both 
massive and tracer) and the integer, m, is the azimuthal 
wave number. Bulge and halo particles are ignored. We use 
logarithmically spaced radial bins. We compute the Fourier 
transform of the complex function composed of Wc"(r, t) 
and W]?(r,t) or 

W m {u, r) = f 2 e iut h(t) [WS*(r, t) + iW^(r, t)] dt (2) 

where h(t) is a Hanning function spanning our time window 
[Ti,T2] between snapshot time T\ and T%. In our spectro- 
gram figures we show the amplitude of the complex Fourier 



components or |W /m (cj,r)|. In our spectrograms the angular 
frequency, uj, increases along the y-axis and log 10 r increases 
along the x-axis). 

A bar or two-armed spiral structure has strong m — 2 
Fourier components in its density distribution (calculated 
with equation [T]) . The pattern speed of the wave is related 
to its angular frequency by £l p = uj/m. Pattern speeds are 
found by dividing the angular frequencies of horizontal fea- 
tures in the spectrograms by the integer m. 

During an N-body simulation with a non-stationary or 
live halo, lopsided modes or waves may develop and the 
bulge and central disk of the galaxy may not remain at a 
fixed position. If the Fourier components are measured re- 
spect to a radial position that is not centered on the bulge 
centroid then an m = 1 component can generate other 
Fourier components. This is sometimes called 'aliasing' or 



'mode-mixing' (e.g., Oiling & Dehnen 20031. As a conse- 
quence we subtract the position of the centroid of the galaxy 
bulge prior to computing the Fourier coefficients with m > 1. 
We will also discuss the m = 1 spectrogram computed with- 
out subtracting the bulge centroid position. 



4.2 Pattern speeds 

After T ~ 400 Myr the bar becomes stable and then evolves 
more slowly during the rest of the simulation. Well defined 
frequencies are seen in spectrograms after bar formation. 
Figures [4] shows spectrograms constructed from m = 2 
Fourier components, whereas Figure[5][6]and[7]show spectro- 
grams constructed from m = 4, 3 and 1 Fourier components, 

is measured 
is measured 



respectively. The spectrogram from Figure 

from the middle of the simulation whereas 

from the second half of the simulation. The same time peri- 
ods are used for Figures[5]-[7| We can see from a comparison 
of Figure and Figure |4Jd (also see Figures [l] and |3| that 
the bar grows longer and its rotation rate slowly decreases 
during the simulation. At 400 Myr the bar is only 6 kpc long 
(semi-major axis) but by 700 Myr its semi-major axis has 
increased to 11 kpc. 

Pattern speeds of waves can be estimated from a spec- 
trogram by dividing the angular frequency of a strong fea- 
ture by the integer m used to make the spectrogram. In our 
spectrograms we plot the angular frequency, u, rather than 
the pattern speed Q = w/m, following most previous works 
(e.g, |Sellwood fc Sparke|1988||Masset fc Tagger|1997[ ). Esti- 



mated pattern speeds for waves seen in the presented spec- 
trograms are listed in Table [I] The strongest pattern seen 
in the m = 2 spectrogram shown in Figure |4^i is the bar 
wave that has an angular frequency of ujb ~ 0.080 Myr -1 
corresponding to a pattern speed of Q.b ~ 0.040 Myr - . 
There are slower patterns at larger radii seen both in the 
middle and later parts of the simulation. These correspond 
to two-armed spiral patterns. Figure [6] showing the m = 3 
spectrograms shows that there is also at least one strong 
three-armed pattern present in the simulation. 

Figure [4] and [5] show curves corresponding to the angu- 
lar frequencies m(f2±«/2), m(Sl±«:/4), mf2 with m — 2 and 
m — 4, respectively where m is that used to compute the 
spectrogram. Each pair of lines (when there is a ± in the 
above listed frequencies) is the same color in the Figures. 
These lines can be used to identify inner and outer Lindblad 
Resonances. Figure [6] and [7] in addition show curves corre- 
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sponding to the angular frequencies m(fi±/t/3), with m = 3 
and m = 1, respectively. Here fi(r) and «(r) are the angu- 
lar rotation rate and the epicyclic frequency. We computed 
these frequencies using a rotation curve generated from the 
azimuthally averaged midplane gravitational potential mea- 
sured from the positions of all the massive particles at time 
T = 500 Myr. 

A Lindblad resonance (LR) with a pattern with speed, 
Q p , and integer m occurs where 



m(fl — Q p ) = ±/t 



(3) 



( |Lindblad|1926t . The sign in the above determines whether 
the resonance is an inner or outer Lindblad resonance (ILR 
or OLR). A corotation resonance (CR) occurs where Q, = 
fi p , or the angular rotation rate equals the pattern speed. 
Here we refer to each resonance by its integer, here denoted 
m. Other works refer to the m — 4 ILR as the 4:1 ILR 
(e.g., |Meidt et al.||2009| ) or the 4/1 ILR (e.g., |Rautiainen| 
& Salo 19991. We also use the notation ILR4 to refer to 
the m = 4 ILR. The radii of the estimated Lindblad and 
corotation resonances for pattern speeds estimated from the 
spectrograms are also shown in Table[l] Resonance locations 
are estimated using the rotation curve computed from the 
azimuthally averaged density distribution at time T = 500 
Myr. The bar pattern also contains non-zero Fourier com- 
ponents in its density distribution with m / 2. As the bar 
lacks symmetry when rotated 180°, or when reflected about 
its major or minor axis, it also contains non-zero odd-m 
Fourier components in its density distribution. 

We see in Figure [4^l that the bar pattern extends out to 
its corotation radius where the bar pattern intersects 2S1_b 
(the dark blue line) at about 6 kpc. The pattern is seen past 
the bar's corotation radius in the spectrogram, consistent 



with previous studies (e.g., see Figure 3 by Rautiainen & 
Salo||1999 L Between the corotation resonance and the bars 
m = 2 outer Lindblad resonance (OLR2), periodic orbits 
can be elongated and oriented perpendicular to the bar. 

Recent studies of the Milky Way have revealed two bar 
like structures, a long thin bar (Benjamin et al. 2005| ) at 
about 45° from the Sun/Galactic center line and half length 
4.4 kpc and a shorter or traditional thicker bar (or triaxial 
bulge component) at about 15° ( Vanhollebeke et al.| |2009). 
Estimates of the bar's pattern sp eed are in the range 50-55 



km s" 1 kpc" 1 or ~ O.OSSMyr" 1 | Dehnen||2000 



Minchev et 



aT1|2007| |Gardner et al.||2010[ ). As shown by |Gardner et al. 
(20101 the Hercules stream is consistent with either long or 



short bar as long as the pattern speed places the solar neigh- 
borhood just outside the bar's OLR2. Our simulated bar has 
with half length ~ 5 kpc and so is longer than the Milky 
Way's long bar but this is consistent with its slower pattern 
speed at 0.040 Myr -1 or 40 km s" 1 kpc -1 . In our simulation 
a position corresponding to that of the solar neighborhood 
would be at a galactocentric radius of about 10 kpc, just 
outside the bar's OLR. 

At both earlier and later times, features at lower pattern 
speeds are also seen in the m = 2 spectrograms. A peak 
just past the end of the bar is seen with a pattern speed of 
fl s ~ 0.030 Myr" 1 at earlier times (Figure^) and at Q s ~ 
0.022 at later times (Figure^). This would correspond to 
an inner two-armed spiral pattern. This pattern appears to 
end within its own corotation radius at 7 or 9 kpc (early and 
late, respectively). There is also a slower outer two-armed 



spiral pattern primarily seen at early times with a pattern 
speed Q s ~ 0.015 (though possibly also at later times with 
tt s ~ 0.010). This pattern also ends near its own corotation 
resonance. 

There is a weak and narrow higher frequency peak in the 
m — 2 spectrograms near the end of the bar at an angular 
rotation frequency of lo ~ 0.13 and co ~ 0.11 Myr" 1 at early 
and late times, respectively. Unlike the slower frequency in 
the m — 2 spectrograms that end near their corotation radii, 
this one extends between its corotation radius and its outer 
Lindblad resonance. There is a similar higher frequency fea- 
ture in the m = 4 spectrograms (Fig. |5| at cj ~ 0.23 and 
0.17 Myr" 1 at early and later times, respectively. 

The m = 4 spectrograms shown in Figure [5] look re- 
markably similar to the m — 2 spectrograms (Figure|4]| . Two 
armed and oval density distributions often contain m — 4 
components so this is not necessarily surprising. The higher 
frequency features near the end of the bar are also seen in 
the m = 4 spectrograms. With the exception of the bar it- 
self, the peaks in the m = 4 spectrogram extend between 
their inner and outer m = 4 Lindblad resonances (shown 
as green lines). This differs from the two-armed spiral pat- 
terns that end within their corotation radii. In the m — 4 
spectrogram the inner spiral pattern appears split into two 
different peaks. We may be able to better measure angular 
frequencies in the m = 4 spectrogram than the m = 2 spec- 
trogram because the angular frequencies are higher and so 
less affected by the precision limit set from the width of the 
time window used to compute the spectrograms. 

Unlike the inner and outer two-armed patterns, the 
three-armed pattern both at early and later times extends 
through its corotation resonance. This pattern (see yellow 
lines in the m = 3 spectrograms shown in Figure [6f extend 
between its m = 3 inner and outer Lindblad resonances. The 
pattern speed is difficult to measure at early times because 
there are two nearby features in the spectrogram. Both at 
early and later times the pattern speed of the three-armed 
spiral is about fi s ~ 0.023 Myr" 1 . 



As noted by Shaviv ( 2003 1 (see their Table 3) pattern 



speeds estimated for spiral patterns in the solar neighbor- 
hood fall into two groups. One group has pattern speed ~ 27 
km s" 1 kpc" 1 similar to the angular rotation rate of a parti- 
cle in a circular orbit, Qq, at the solar galactocentric radius, 
Rq, and the other group has a pattern speed ~ 14 km s" 1 
kpc" 1 . The mock solar neighborhood position in our simula- 
tion lies outside the inner two armed spiral pattern at 0.030 
Myr" 1 , but is affected by the three armed spiral wave with 
pattern speed 0.023 Myr" 1 or ~ 23 km s" 1 kpc" 1 and the 
outer two armed spiral with a pattern speed of 0.015 Myr" 1 
or ~ 15 km s" 1 kpc" 1 . If the spiral waves present in the solar 
neighborhood are driven by the Galactic bar then it might 
be interesting to match the faster pattern speed measured 
in the solar neighborhood with a three armed pattern driven 
by the bar and the slower pattern speed with a two armed 
pattern also possibly linked to the bar. However, with the 
exception of |Naoz fc Shaviv] ( |2007| ), few studies have con- 
sidered interpretation of the solar neighborhood in terms of 
two or more patterns. 

Previous works have seen spiral structure in N-body or 
particle mesh simulations at pattern speeds that differ from 
the bar's (e . g., |Sellwood fc Sparke||1988| [Patsis fc . Kauf-| 
mann||1999| |Rautiainen fc Salo||1999| |Voglis et al.||2006[) 
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Masset & Tagger (19971 interpreted this behavior in terms 
of non-linear mode coupling ( Tagger et al.f l987 Sygnet et 
|aT1|1988| ), whereas |Sellwood & Lin|(|1989|) interpret the dif- 



with m = 4 ( jMasset fc Tagger|1997| . Like |Masset fc Tagger 
( 1997 1 we see features in the m = 4 spectrograms at fre- 



ferences in pattern speeds in terms of transient behavior and 
a groove instability. The spiral patterns in the simulations 
by Sellwood & Sparke ( 1988 1 extend from their m = 2 in- 



ner Lindblad resonances to their outer Lindblad resonances 
(see their Figure 1). However, the patterns of the two spiral 



density waves seen by |Patsis fc Kaufmann ( 1999 1 in a simu- 
lation lacking a bar extend from their m = 2 inner to their 
outer m = 4 Lindblad resonances. Both of these situations 
differ from the two-armed structures in our simulations that 
end within their corotation radii. The simulations by|RautL"| 



ainen & Salo ( 1999 1 are closest in behavior to ours, showing 



multiple spiral patterns extending past the bar, with some 
patterns ending near or within their corotation radii. 

As the spiral and bar patterns in our simulations are 
moderately long lived, the spiral patterns could be driven by 
the bar (|Yuan fc Kuo|p97l |Masset fc Tagger ||l997l |Rauti 



ainen fc Salo| 1999 1. Rautiainen & Salo ( 19991 accounted for 
the ranges and speeds of spiral patterns using the non-linear 
wave coupling scenario. They found that they could often as- 
sociate the beginning of a slower outer pattern with the loca- 
tion of one of its resonances that is also a resonance with an 
inner faster bar or other spiral pattern. The overlap between 
the resonances relates one pattern speed to the other and 
couples one perturbation to the other. jRautiaincn fc Salo| 
(1999) found examples of CR/ILR 4 , CR/ILR2, OLR2/CR 
and OLR2 /ILR2 couplings between spiral and bar waves and 
between two different spiral density waves. 

For our system there is a possible connection between 
the bar's CR and the m = 3 ILR of the three-armed spiral 
armed wave. A comparison of resonance locations (see Ta- 
ble [l| shows that the bar's corotation resonance lies close to 
the three-armed spiral's m = 3 ILR. The radial distance be- 
tween these resonances is small compared to the error from 
measurement in the spectrogram. There is also a possible 
coupling between other resonances. It is difficult to differ- 
entiate between possibilities by resonance locations alone as 
driving could occur in a spatially broad region near a reso- 



nance. As suggested by Masset & Tagger ( 19971 we can test 



for mode or wave coupling by searching for waves at beat 



frequencies though Rautiainen & Salo ( 1999) found that ex- 



pected beat waves are not always detectable. 



4.3 Non-linear Coupling of Waves 



The scenario introduced by Tagger et al. ( 1987 1 suggests 



that there are strong non-linearities in the stellar response 
when resonances coincide. Two waves (e.g., bar and spiral) 
couple to a third one at a beat frequency. For example we 
consider azimuthal wavenumber tub associated with a bar 
perturbation at pattern speed £Ib and so frequency ujb — 
nriB^lB- We consider azimuthal wavenumber ms associated 
with a spiral pattern at pattern speed Q s and so angular 
frequency ujs = ms^s- A third wave can be excited with an 
azimuthal wavenumber m = ms+ms and angular frequency 
uj = ujb + ujs or with wavenumber m — tub — ms and 
uj = lob — uss- 



An example explored by Tagger et al. ( 1987 1 is a bar 



and spiral mode, both with m = 2, coupled to a slow ax- 
isymmetric mode with m — and a faster four-armed one 



quencies above that of the bar, and these can be considered 
predictions of the wave-coupling scenario. Unfortunately the 
number of features present in the m — 2 and m — 4 spectro- 
grams is high and the imprecision of measured frequencies 
makes it difficult to definitively determine whether there is 
coupling between m — 2,4 waves and a slow axisymmetric 
mode. 

As our simulation contains a bar (m = 2), a lopsided 
m — 1 motion and a three armed wave we can consider cou- 
pling between these three waves. To create the m = 2, 3 
and 4 spectrogram we first subtracted the bulge centroid, 
however the galaxy center moves during the simulation. In 
Figure [7j showing the m = 1 spectrograms but constructed 
in the inertial frame, we see features associated with the 
major patterns identified in the other spectrograms. For ex- 
ample we see a peak at uj — 0.08 Myr^ 1 that coincides with 
the bar's angular frequency. One interesting peak is at an 
angular rotation rate of uj ~ 0.010 that is also seen in the 
m = 3 spectrogram at a radius of about 8 kpc (logior = 0.9). 
At earlier times there is power at cj ~ 0.03 in the m — 1 
spectrogram that is broader than the similar feature seen in 
the m = 2 spectrogram. This suggests that there is sufficient 
power at low frequencies in the lopsided wave that it could 
couple to other waves. 

Since there is power at low frequencies in lopsided mo- 
tion, we can consider the possibility of a coupling between 
a slow lopsided mode and the bar. Here the index and pat- 
tern speed of the lopsided wave is rriL = and that of 
the bar is tub = 2, fie. The sum leads to a three-armed 
spiral pattern with angular frequency ujs = £Il + 2CIb or 
pattern speed Qs = § (fii + 2f2_g). If we assume that this 
corresponds to the three-armed pattern with Sis ~ 0.028 
then we can solve for the angular frequency of the lopsided 
mode Q.l ~ 0.01. There may be power at this frequency as 
there are low frequency features in the m = 1 spectrogram. 
Thus, a lopsided slow mode could be coupling the bar and 
the three armed spiral waves. 

A variety of resonant overlaps or non-linear wave cou- 
plings are possible and may account for the driving of the 
patterns present in the simulation. If we could more precisely 
measure frequencies of the patterns then we would better 
determine which associations are most important or likely. 
Unfortunately peaks in the spectrograms are sensitive to 
the range (in time) of simulation data used to construct the 
spectrogram (e.g., compare Figure^ to Figure^). Longer 
time series are required to make more accurate frequency 
measurements. Since the pattern speeds are not fixed, in- 
creasing the size of the time window does not increase the 
precision of a frequency measurement. This is particularly 
the case for slow patterns that might be important for the 
non-linear coupling models. 

In summary, we find that the spectrograms are rich in 
small features. The shape of the spectrograms seen at early 
times resembles that at later times implying that the fre- 
quencies of these smaller features are slowing along with the 
bar pattern. Even though there are many waves present, 
they are not stochastically appearing and disappearing but 
slowly evolving. Lindblad resonances associated with one 
wave are likely to be coincident with resonances of another 
wave present in the simulation. Together these findings sug- 
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gest that the waves present in the simulation are coupled 
due to non-linear interactions. The most convincing and 
strongest coupling is that between the bar, the lopsided wave 
and the three-armed spiral wave. The multitude of small 
features present in the spectrograms suggest that additional 
couplings are present. The presence of features in the spec- 
trograms at frequencies above that of the bar is consistent 
with the wave-coupling scenario. 



5 STRUCTURE IN THE UV PLANE 

We first describe how we construct velocity distributions 
in local neighborhoods. We then discuss possible relations 
between features seen in these velocity distribution and the 
spiral or bar patterns. 



5.1 Velocity uv distribution histograms in local 
neighborhoods 

We construct velocity distribution histograms in different 
local neighborhoods in the following manner. Consider a 
neighborhood centered about a point in the disk plane 
with coordinates (20,2/0,0) or in cylindrical coordinates 
(ro,</>o,0). We will refer to neighborhoods by the cylindri- 
cal coordinates of their centers; (ro, 4>o)- We consider a star 
or simulated particle in the neighborhood of this center 
point if its coordinates (x, y, z) have distance in the plane, 
R — y/ (x — xo) 2 + (y — yo) 2 < ro/ where / is a dimension- 
less scaling factor setting the size of the neighborhood. We 
compute the radial and tangential components of the veloc- 
ity of each star with respect to the vector to the galactic 
center defined by the central coordinates of the neighbor- 
hood, (a;o,J/o,0). The u component is —1 times the radial 
velocity component in this galactic coordinate system. The 
v component is the tangential component in this coordi- 
nate system minus the circular velocity estimated from the 
azimuthally averaged midplane rotation curve computed at 
time T = 500 Myr. 

Histograms showing the uv velocity distribution are 
computed from the disk particles in each neighborhood. The 
computed histograms are effectively coarse-grained phase 
space distributions. The larger the neighborhood, the larger 
number of stars in the neighborhood. If there are too few 
stars in the neighborhood then the velocity distribution is 
noisy. If the neighborhood is large then structure present 
on small spatial scales may be masked or smoothed by the 
coarse graining. An extreme example would be if the neigh- 
borhood were so large that it contained two spiral arm den- 
sity peaks. 

A series of velocity distributions (uv histograms) is 
shown in Figure [8] at simulation time T =600 Myr. Velocity 
distributions are extracted from neighborhoods centered at 
8 equidistant angles and at 6 radii spaced logarithmically. 
Each row contains histograms from neighborhoods at a par- 
ticular galactocentric radius. The lowest row has ro = 4.1 
and the top row ro = 12.5. Each radius differs from the pre- 
vious by a factor of 1.25. Each column shows neighborhoods 
with a particular galactocentric angle with the leftmost col- 
umn aligned with the bar. Each neighborhood is separated 
in angle by 45° with angle increasing counter-clockwise in 
the galaxy plane. Each neighborhood has a size set by the 



factor / setting the maximum distance of stars included in 
the neighborhood from its center. For the neighborhoods 
shown in Figure [8] the ratio of the neighborhood radius to 
ro is / = 0.1. Positions and sizes of each neighborhood are 
shown in Figure [9] in both Cartesian and polar coordinates. 
At ro = 8 kpc the radius of the neighborhood corresponds 
to 0.8 kpc which is small compared to some N-body studies 
(e.g., Gomez et al. 2010 who used a neighborhood radius 



of 2.5 kpc) but larger than the distance of stars in solar 
neighborhood samples (e.g., Siebert et al.|2011 |. 

We have compared the velocity distributions created 
with / = 0.1 and with / = 0.05 with smaller neighbor- 
hoods. A comparison at time 800 Myr and at a radius of 10 
kpc is shown in Figure [To| We also show a velocity distribu- 
tion constructed of massive disk particles only to illustrate 
that we see no evidence for significant differences between 
massive and massless disk particle populations. The velocity 
distributions in smaller neighborhoods are noisier as would 
be expected since each neighborhood contains fewer parti- 
cles. The clumps in the plots are slightly sharper, otherwise 
the morphology is almost identical. The density distribution 
is relatively smooth (see Figure[9| implying that the velocity 
dispersion is not low in the disk (though it is lower at larger 
radii than at smaller radii). A simulation of a colder disk 
with more tightly wound spiral arms might show features in 
the velocity distributions that are more strongly dependent 
on the size of the neighborhood. 

In Figure [8] we see that multiple clumps exist in the lo- 
cal velocity distribution at all radii in the galaxy. Velocity 
clumps were previously seen in a local neighborhood of the 



N-body simulation presented by Fux ( 2000 1. Since our simu- 



lation lacks mergers, the clumps in the velocity distribution 
are caused purely by dynamical phenomena in the disk (as 
also concluded by Fux 20001. Clumps in the velocity dis- 



tribution in one neighborhood are often also present in the 
nearby neighborhood at a larger or smaller radius but at a 
shifted v velocity. For example at an angle of 315° there is an 
arc with a mean velocity of v ~ 50 km/s for neighborhood 
radius ro = 4.1 that has a lower velocity of v ~ 10 km/s 
for at ro = 5.1. We can consider the possibility that a single 
halo particle could be responsible for one of these features. 
The separation between the two neighborhoods mentioned 
above is a kpc. However the circular velocity at a distance of 
a kpc from one of our halo particles is only 2 km/s. Thus a 
single halo particle passing through the disk cannot account 
for the correlated and broad structures seen in the veloc- 
ity distributions. We will discuss the relationship between v 
and the mean orbital or guiding radius below when we con- 
sider the relation between velocity components and orbital 
properties. 

5.2 Velocity distributions in arm and interarm 
regions 

When looking at the velocity distribution in different neigh- 
borhoods it is helpful to know which neighborhoods have a 
high surface density. Figure [9^ and b also show the loca- 
tions of neighborhoods on the projected disk surface density 
used to construct the velocity histograms shown in Figure 
[8] In Figure , projected in polar coordinates, the x axis 
gives the azimuthal angle with respect to the bar. Thus each 
neighborhood shown on this plot directly corresponds to the 
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position of the uv histogram in Figure [8] This makes it pos- 
sible to see how the velocity distribution varies with position 
and between arm and interarm regions. 

A comparison between Figure [8] and [9)3 allows us to 
see if the velocity distribution in high density regions dif- 
fers from that at low density regions. We find that when 
the density increases the velocity dispersion is also likely to 
increase. The velocity distribution is wider when the neigh- 
borhood lies on an arm or bar peak than in interarm regions. 
Phase space density should be conserved as particles orbit 
the galaxy, so this trend is expected. 

We can compare the velocity distribution on the leading 
(or pre-shock) side of a spiral arm to that on the trailing 
(or post-shock) side. Our simulated galaxy rotates counter- 
clockwise or with increasing (j>. We denote the leading side 
of a spiral arm as that located at larger <f> than the peak, 
or to the right in Figure [9}d (and counter clockwise in angle 
from an arm peak in Figure [9^). Observed galaxies tend 
to have dust lanes on the trailing side (as denoted here) 
of a spiral arm. This happens because the spiral pattern is 
usually faster than local galactic rotation (in other words 
the corotation radius lies interior to the pattern) . In a frame 
moving with the gas on Figure [9Jd, the pattern moves to the 
right and so shocks would form on the left-hand or trailing 
side of each density peak. 

We see that arcs in the velocity distribution with posi- 
tive u and v are often seen on spiral arm density peaks and 
on the leading (or pre-shock) side of an arm. These arcs are 
oriented with negative slope such that v increases with de- 
creasing u. Moving from the leading to the trailing side at 
a particular radius (comparing histograms from right to left 
in Figure [sj the high v arcs move from positive to negative 
u. See for example r = 8 and (j> — 90, 135° or r = 6.4 and 
<j) = 315,0° in Figure [8] Interarm regions tend to lack arcs 
that are centered at high and positive u, v. 

In interarm regions the velocity dispersion is lower and 
there is sometimes a clump, separated from the dominant 
peak, that is at negative u and v. The clump can be elon- 
gated in u, resembling the Hercules stream in the solar 



neighborhood's velocity distribution (e.g., Dehnen & Bin- 



ney| 1998| [Arifyanto fc Fuchs|2006||Gardner et al.|2010||Bo~ 
20101). 



5.3 Gaps in the uv histograms 

Figure [12] shows the neighborhoods at a radius of ro = 6.4 
kpc for the same angles shown in Figure [8] but for snapshots 
separated by 50 Myr between 500 and 950 Myrs. The angles 
here are given with respect to the bar and so in a coordinate 
system rotating with the bar. At each angle with respect to 
the bar, clumps do not remain fixed as the system evolves. 
This implies that spiral structure influences the locations 
of gaps and clumps in the velocity distribution. The gaps 
cannot be solely due to bar perturbations. 

A stream similar to the Hercules stream at v ~ — 50 
km/s is present in the velocity distribution at early times 
at angle (aligned with the bar) at ro = 6.4 kpc. Like 
the Hercules stream (e.g., [Deh nen 1999} this arc or clump 
has a wide u distribution, a narrow v distribution and is 
centered at a negative u value. However the bar's m = 2 
OLR (at about 9 kpc) lies outside these neighborhoods, so 
this clump cannot be associated with this resonance, Similar 



feature low u, v clumps in the velocity distribution are also 
seen at other radii such as ro = 8 and 10 kpc but appearing 
at different times in the simulation and at different angles 
with respect to the bar, consequently the bar and resonances 
associated with it cannot provide the only explanation for 
these features. 

The spiral density waves have slower patterns than the 
bar, hence density peaks associated with them move to the 
left (or to a smaller <j>) as a function of time or moving 
downward on Figure |12[ At one position as a wave peak 
passes through the neighborhood first a velocity distribution 
typical of a leading side of an arm is seen then one typical of 
a trailing side. The velocity distribution contains a positive 
u and v feature when on the leading side of a spiral feature. 
The feature is centered at progressively lower u values as the 
arm passes by and the neighborhood moves to the trailing 
side. See for example the right most column and at an angle 
of 315° between times 500 to 700 Myrs (top 5 rows). 

While clumps come and go in time, we find that the 
gaps between the clumps tend to remain at the same v val- 
ues. For example consider the right most column of Figure 
|12[ Two gaps are seen with different degrees of prominence 
at differ times. We recall from Figure [8] that gaps at one 
neighborhood radius were often seen at another neighbor- 
hood radius but shifted in d. To compare a gap seen in a 
neighborhood at one radius to that seen in a neighborhood 
at another radius we can consider motions of particle orbits. 



6 INTERPRETING STRUCTURE IN THE UV 
PLANE 

To relate structure in local velocity distributions it is help- 
ful to recall models for orbits. We first consider epicyclic 
motion in the absence of perturbations, then review the or- 
bital dynamics to first order in the perturbation strength in 
the presence of perturbations. In both case we consider the 
relation between uv velocities in a local neighborhood and 
quantities describing the orbit such as a mean or guiding 
radius and an epicyclic amplitude. 

In the absence of perturbations from spiral arms, the 
motion of stars in the disk of a galaxy can be described 
in terms of radial or epicyclic oscillations about a circular 
orbit ( |Lindblad|[T926[ |Kalnajs||1979[ ). It is useful to specify 
the relation between the observed velocity components u, v 
and the parameters describing the epicyclic motion (e.g 



Fux 



[2000} [200i] |Quillen fc Minchev|200"5] > or the mean radius or 
guiding radius, r g , and the epicyclic amplitude. The energy 
of an orbit in the plane of an axisymmetric system (neglect- 
ing perturbations from spiral structure) with a flat rotation 



E(u,v) 



(v c +v) 2 u 



+ — + v c In r + constant 



(4) 



where the potential energy, lnr, is that appropriate for a fiat 
rotation curve, r is the Galactocentric radius and v c is the 
circular velocity. 

With an epicyclic approximation we can write the en- 
ergy 



E = vi 



+ E e 



(5) 



where the term on the left is the energy of a star in a circular 
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orbit about a guiding radius r g and E ep i is the energy from 
the epicyclic motion, 



2/ \2 2 2 

K \r-r g ) _ K a , > 

2 ' 2 (6) 

Here a is the epicyclic amplitude and k is the epicyclic fre- 
quency at the guiding radius r g . 

We now consider stars specifically in a neighborhood 
of galactocentric radius r , restricting us to a specific loca- 
tion in the galaxy. Setting the energy equal to that written 
in terms of the epicyclic motion using equations (15) [m , we 



solve for the guiding radius, r g . It is convenient to define 
the distance between the guiding or mean radius and the 
neighborhood's galactocentric radius, s = r g — ro. To first 
order in v and s we find that a star in the neighborhood 
with velocity components u, v has a guiding radius with 



ro 



and epicyclic amplitude a 



ro 



U , 2 



(7) 



(8) 



The factor of a half is from k 2 /O. 2 equivalent to 2 when the 
rotation curve is flat. An angle describing the phase in the 
epicycle is 



at an 



V2i 



(9) 



(consistent with equations 2-4, Quillen|2003 for a flat rota- 
tion curve). The epicyclic phase angle = at apocenter. 

The above three relations allow us to relate the veloc- 
ity components it, v for stars in the solar neighborhood, to 
quantities used to describe the epicyclic motion; the guiding 
radius, epicyclic amplitude and phase angle. Particles with 
positive v have s > and so guiding radii that are larger 
than the radius of the neighborhood, ro. Particles with neg- 
ative v have s < and so are expected to spend most of their 
orbits inside ro. The distance from the origin, u = v = 0, 
determines the epicyclic amplitude. 

It may also be useful to write 



JL 
ro 



(10) 



where we have used our definition for s in equation [7] This 
equation shows that to first order v sets the guiding radius 
of a particle. This is useful because the guiding radius sets 
the approximate angular rotation rate and rotation period 
of the orbit. The guiding radius also sets the approximate 
epicyclic frequency (or period) of the particle's orbit. Reso- 
nances occur where sums of integer multiples of these peri- 
ods are equal to an integer multiple of a pattern speed. Thus 
the v value of clumps or gaps in the velocity distribution may 
also allow us to locate resonances. 



6.1 Spiral Perturbations to First Order in the uv 
plane 

The above approximation neglects the effect of bar or spiral 
perturbations. Spiral waves give a time dependent perturba- 
tion to the gravitational potential that can be approximated 
by a single Fourier component depending on a single fre- 
quency. For a logarithmic spiral perturbation with m arms 



we can consider V(r, 9, t) = V m cos(m(9—Q p t+a In r)) where 
the pitch angle a = —d8/dr. To first order in the perturba- 
tion strength V m the orbit of a star in the galactic midplane 
with mean or guiding radius r g 

r(r g ,9,t) = r g + acos(K g t + cj>o) + 

C(r g ) cos(m(9 — Q p t + alnr 3 + cq)) 



(11) 



( Binney fc Tremaine||1987 equation 3- 119a) where a is an 
epicyclic amplitude, and Co, </>o are phase angles. The forced 
epicyclic amplitude in the WKB approximation |C(r 9 )| « 

where A s — k 2 — m 2 (Q. g — Q p ) 2 represents the dis- 
tance to a Lindblad resonance and is evaluated at the guid- 
ing or mean radius, r g and k 9 and Q g are the epicyclic fre- 
quency and angular rotation rate at r g . The denominator 
A g becomes small near the Lindblad resonance and the first 
order approximation breaks down. A higher order approxi- 
mation can be used to show that closed orbits exist on both 
sides of resonance and that epicyclic amplitudes do not be- 
come infinite on resonance (e.g., Contopoulos|1975 1. 

The orbital motion (equation 11 1 consists of an epicyclic 
motion with amplitude a and a phase <j>o and a forced motion 
aligned with the potential perturbation and moving with it. 
A population of stars would have a distribution in both a 
and 0o. When a — and m — 2 the orbits are ellipses. The 
orbit with a = is periodic or closed in the frame rotating 
with the spiral pattern. The mean value, a, sets the velocity 
dispersion in the disk. When a < \C\ the orbits are nearly 
closed in the frame rotating with the potential perturbation 
and we can consider the disk to be cold. 

The radial and tangential velocity components can be 
computed 



v r (r g ,l 
Ve(r g , I 



.*) 
.*) 



—CbK g sm(Kgt + cj>o) (12) 
—C(r g )m(£l g — Q p ) sin(m(# — Q p t + a\nr g + Co)) 
r g Sl g — aQ g cos(frf + <j>o) (13) 
— C(r g )Q g cos(m(9 — Q. p t + a\nr g + Co)). 

(from equation 3.117 by Binney & Tremainc 1987]) . When 
the rotation curve is flat r g Q g is independent of radius and 
the uv velocity components u = — v r and v — vg — r g Q g . 

We note that the equations for r(r g ,9,t) and ve(r g ,9,t) 
contain similar terms but with the opposite sign. A shift in v 
corresponds to a shift in radius and the mean orbital radius 
can be roughly estimated from the v component as was true 



in the case of pure epicyclic motion (equation 10 1. If the free 
epicycle is small (a = 0) or the disk is cold then 



U l ilg 

v \m(Q — fl p ) 
with 



- tan md 



6i = (9 — Q p t + a In r g + c ) 



(14) 



(15) 



and is sensitive to the angle with respect to the potential 
perturbation (compare to equation[9] for pure epicyclic per- 
turbations). Likewise the distance from origin on the u,v 
plane is dependent on the amplitude of radial oscillations 
with 



u 2 Q, 2 



n p y 



C ^ BJ " V m 2 (ft g 
(compare to equation [8] 



+ v 2 



(16) 
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In the limit a < \C\, particles are in nearly elliptical 
orbits but the orientation of the ellipses can vary with ra- 
dius. A set of ellipses with differing orientations is shown 
in Figure [TT^ . Each ellipse has the same eccentricity but is 
shifted in angle with respect to the previous one. The differ- 
ence in angle is 20° except between the 6th and 7th ellipse. 
The difference in orientation angles between these two el- 
lipses is larger, 90°, and introduces a discontinuity that we 
will discuss later when we consider the relationship between 
gaps in the velocity distribution and discontinuities in the 
spiral structure. The red circle on Figure [TT^ , shows an ex- 
ample of an arm peak that lies outside the discontinuity. In 
this neighborhood there are nearly overlapping orbits. The 
range of velocities seen in this neighborhood would be high 
as there are stars from many different orbits passing through 
the neighborhood. There should be a relation between u and 
v velocity components as there is a relation between mean 
galactic radius of the ellipse and orbital angle in the neigh- 
borhood. Thus the velocity distribution would exhibit an 
arc. 

The green circle in Figure [TT^ shows an interarm region 
where the angles of orbits only varies slightly compared to 
that in the red neighborhood. In the green interarm neigh- 
borhood we expect a narrow velocity distribution. This il- 
lustrates why we see expect to see arcs in the velocity distri- 
butions on spiral arm peaks and lower velocity dispersions 
in interarm regions. 

An arc in a local velocity distribution that has u de- 
creasing as ii increases can be interpreted in terms of orbital 
properties as a function of r g (approximately set by v). The 
positive v orbits have larger guiding radii. The change in 
angle or v ju slope on the velocity distribution plots implies 
that the epicyclic angle <f>i (equation 15 1 decreases with in- 
creasing r g . Thus these arcs correspond to a relation be- 
tween epicyclic angle and guiding radius. Consider our car- 
toon, Figure [TT^,, showing concentric ellipses with different 
orientation angles. A neighborhood on top of a spiral arm 
peak (such as the red one) intersects orbits with a range of 
epicyclic angles, corresponding on this figure to a smooth 
change in orientation of the ellipse. Thus we expect that an 
arc would be observed in a velocity distribution at the posi- 
tion of a spiral arm peak. The angle of the arc in the velocity 
distribution is likely to be affected by the winding angle of 
the spiral structure with a more open spiral corresponding to 
a arc with steeper slope. However a neighborhood may con- 
tain orbits with a smaller range of guiding radii if the spiral 
structure is less tightly wound. Model orbital distributions 
would be needed to use the slopes and ranges of velocities 
of clumps seen in velocity distributions, and their gradients, 
to place constraints on the spiral structure (pattern speed, 
pitch angle and strength) from observations of the velocity 
distribution alone. 

We now consider the type of arc that would be seen 
on a trailing spiral arm. The red neighborhood in Figure 
|f I Y intersects orbits near apocenter with small guiding radii 
(and so negative v) and orbits that are near pericenter that 
have larger guiding radii (and so positive v). In between and 
in the neighborhood are orbits that are moving toward the 
Galactic center and so have positive u values. For trailing 
spiral structure we expect the arc to progress from low v 
to high v passing through the top right quadrant of the uv 
plane, as we seen in arm peaks in Figures [8| |I2| and |I3| Thus 



our illustration Figure |II^ provides an explanation consis- 
tent with spiral structure for the arcs seen in the velocity 
distribution. 

We can use our first order perturbation model (equa- 
tions If - 13 1 to consider particles with nearby guiding radii. 



We expand equation (III for closed orbits (a — 0) 



r(r g +x,6,t) w r g + x (17) 
+C(r g ) cos(m(# — fipi + a\nr g + Co)) 
—C(r g ) sin(m(# — Sl p t + a lnr g + Co)) 



where we have used the WKB approximation by neglecting 
dC(r g )/dr. The derivative of the above equation 



dr 
dr„ 



C(r g ) sin(ra(# — Q, p t + a\ar g + Co))- 



(18) 



Consider two closed orbits separated by a small differ- 
ence in guiding radius. The radial distance between these 
closed orbits at a particle angular location is smallest when 
the above derivative is a minimum. A small difference in cur- 
rent radius is possible when the sine term is negative. This 
happens midway between pericenter and apocenter. Thus 
the closed orbits are nearest to each other when this sine 
term is negative (mcj>\ — 7r/2). The radial velocity is also 
proportional to the sine of the angle. For a > 0, correspond- 
ing to trailing arms, we find that the orbits are closest to- 
gether when the radial velocity component is negative and 
the particles are approaching pericenter (positive u, zero v), 
as shown by the red circle in Figure |ri> . Orbits are maxi- 
mally distant when the radial velocity is positive or 90 de- 
grees away from the red circle on Figure 1 1 l^ t, that shows 
a two-armed structure. At locations and corresponding ve- 
locities where the above derivative is high, the phase space 
density would be higher and so there would be features seen 
in a local velocity distribution. Inclusion of the the deriva- 
tive dC(r g )/dr when estimating the derivative dr/dr g would 
shift the estimated angle of the minimum of the derivative. 
This shift may explain why arcs in the velocity distribu- 
tion are primarily located in the upper right hand quadrant 
(positive u, v) and not centered at it > 0, v = 0. 

The above discussion has neglected the dispersion of 
orbits about the closed or periodic ones shown in our car- 
toons (Figure 1 1 1 [ ), however if the dispersion is neglected 
then the predicted velocity distribution contains particles at 
only a few velocities. The stellar dispersion has been used in 
constructing weights for orbits to populate the velocity his- 
tograms and predict velocity distributions from test particle 
integrations (e.g., |Dehn"en]|2000| [Quillen fc Minchev||2005[ |. 
We should remember that Figure [Tl| only qualitatively gives 
us a relation between orbital and velocity distribution, how- 
ever the properties of closed orbits are relevant for interpre- 
tation of these distributions (e.g, Quil len fc Minchev|2005 1. 

The above equations are given in the WKB approxima- 
tion for a spiral perturbation. For a bar or two-armed open 
spiral perturbation to first order the orbits are also ellipses 
however the WKB approximation cannot be used. For a bar 
perturbation an additional term in ve should be taken into 
account when estimating the guiding radius and the angle 
m(6 — Q p t) from u, v (see 



Binney & Tremaine 1987 



section 

3.3). For a spiral perturbation if the WKB approximation 
is not used additional terms in both velocity components 
would be used to estimate the guiding radius, the and orbit 
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orientation and the angle with respect to it from u,v (e.g., 
Contopoulos|1975 l. 



7 GAPS AT DIFFERENT RADII 

We now look at the location of gaps in the velocity distribu- 
tion for the snapshot shown in Figures [8] and [9] at T = 600 
Myr. In Figure [8] at a radius of ro — 5.1 kpc and an angle of 
4>o = 135° there are two gaps seen in the velocity distribu- 
tion, one centered at v ~ 50 and the other at v ~ km/s. 
The one at v ~ has a guiding radius of r 9 ~ 5.1 kpc, the 
same as the neighborhood radius. Using equation |10[ we can 
predict the v value for this gap were it to be observed from 
at a neighborhood radius of ro = 4.1 kpc. We estimate that 
it should lie at v ~ 50 km/s. There is an arc with a gap at 
about v ~ 50 km/s in the velocity distribution at ro = 4.1 
kpc and at 4>o = 135°, as expected. The other gap in the 
neighborhood with ro = 5.1, <f>o = 135° at v ~ 50 km/s has 
a guiding radius r 9 ~ 6.3 kpc (again using equation 10 1. A 
clump with this guiding radius would lie at v ~ km/s in 
a neighborhood with ro = 6.4 kpc. We don't see a gap with 
this v value at ro = 6.4 kpc and 4>o = 135° but we do see one 
in the neighboring position at ro = 6.4 kpc, <j>o = 90°. Thus 
we find that relation between v and guiding radius is con- 
sistent with shifts in individual features in v as a function 
of neighborhood radius. 

We also find gaps at nearly zero v values in neighbor- 
hoods with ro ~ 8 and 10 kpc. In the inner galaxy there is 
a gap at ro = 4.1 kpc, cfio = 45, 90° with slightly negative v 
that is also present at v ~ — 50 km/s in the neighborhood 
with ro = 5.1 kpc, <j>o = 45° and has r 9 ~ 3.9 kpc. Alto- 
gether we see gaps at the following guiding radii r g ~ 3.9, 
5.1, 6.3, 8 and 10 kpc at this time in the simulation. 

Looking back at Figure [8] we compare the location of 
gaps in the velocity distribution (identified by their guiding 
radii) and the locations of kinks (changes in winding an- 
gle) or discontinuities in the spiral arms. In Figure [9J) we 
have identified radii at which we find discontinuities in the 
spiral structure, either changes in winding angle (kinks) or 
places where spiral arms or armlets appear or disappear. 
These radii are shown as black circles in Figure [9^, and as 
horizontal lines in figure [9Jd. The radii of the discontinuities 
on these plots are at rd ~ 3.5, 5.1, 6.3, 8.7 and 11 kpc. A 
comparison between the list of discontinuities and our list 
of gaps identified in the velocity distributions (identified by 
their guiding radii) suggests that these two lists are related. 

A gap in the velocity distribution occurs in a neighbor- 
hood when there are two groups of stars each with velocity 
vectors in similar but different directions, but few stars with 
velocity vectors in between those of the two groups. In Fig- 
ure |llfc , we show concentric ellipses that increase in radius 
and shift in orientation angle. The difference in orientation 
angle between each ellipse is small and constant except be- 
tween the 6-th and 7-th ellipse where the difference is much 
large and 90°. There is a discontinuity in the position of 
the spiral arm peak at the radius of the introduced large 
change in orientation angle that we can see from the shape 
of the orbital overlap regions. The discontinuity introduces 
regions such as the neighborhood illustrated with the blue 
circle. This neighborhood contains orbits with two separate 
groups of orbit angles. In this neighborhood we would ex- 



pect a corresponding gap in the velocity distribution. Figure 
|llfc shows that we expect gaps in the velocity distribution 
where there are discontinuities in the spiral structure, and 
suggests that gaps in the velocity distribution only occur 
at radii where there are discontinuities in spiral structure 
and at specific angles where there are spiral arm peaks. We 
don't see gaps in every velocity distribution in Figure [8] but 
do tend to see them where the surface density is high. This is 
expected as density peaks can occur when multiple pattern 
peaks coincide. 

In the previous section we introduced a first order per- 
turbative orbital model containing a single potential pertur- 
bation (equations 11 - 131. This model can be extended for 



multiple patterns giving orbits with radius, radial and tan- 
gential velocity as a function of guiding or mean radius to 
first order in the perturbation strengths 

r(r g ,9,t) = r g + acos(K g t + cj>o) (19) 
+Ci(r g ) cos(mi(# — flit + ai lnr 9 + ci) 
+C 2 (r 9 ) cos(m 2 (S — fi 2 i + a 2 lnr 9 + c 2 ) 

When a is small compared to |Ci| and | Ca | the orbits are 
nearly closed as shown in Figure |TTfc ,,b. 

For a particular neighborhood with r, 9 the above equa- 
tions must be inverted to solve for r g , a and <po as a function 
of the velocity components. The inversion need not give a 
unique value of r g so particles from separated radial regions 
(in guiding radius r 9 ) can simultaneously be present in one 
local neighborhood. The similarity of terms in r, it, ve im- 
ply that there is a close relation between the inverted r g 
distribution and the velocity components. This situation is 
illustrated in Figures |llfc ,b where transitions between one 
pattern and another gives neighborhoods containing orbits 
from different regions of the galaxy and a bifurcated veloc- 
ity distribution in that neighborhood. Figure 1 1 I^ l shows a 
transition between one two armed and another two-armed 
pattern. Figure |llb illustrates a transition between a two- 
armed and a three-armed pattern. In this figure the inner 
6 curves are ellipses however the outer ones are a sum of 
elliptical and triangle perturbations such as might occur if 
there are simultaneously two-armed and three-armed spiral 
density waves. The elliptical and triangular perturbations 
vary smoothly in orientation angle from curve to curve. The 
transition region exhibits discontinuities in overlap regions 
caused by the onset of the three-armed pattern. The blue 
circle is located at a radius where there is a discontinuity 
and shows that the angular distribution of orbits is bifur- 
cated. We expect that in such a region there will be gap in 
the velocity distribution. 

A neighborhood with r, 6 is most likely to contain two 
separate solutions for r g if the phases of the arguments for 
each perturbation have opposite signs. For example if the 
angle m\(6 — flit + ailnr 9 + ci) ~ and Ci > then 
r > r 9 and the particles in a neighborhood at r will have 
guiding radii interior to r. The orbits will have it ~ but v 
greater than the circular velocity at r (consistent with the 
discussion in section 5.5). If the angle m 2 (# — f2 2 i + a 2 lnr 9 + 
c 2 ) ~ 7r then the opposite is true. A neighborhood at r can 
simultaneously contain stars on orbits associated with an 
interior perturbation near apocenter and stars associated 
with an outer perturbation near peri-center. 

Our spectrograms showed that there are both two- 
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armed and three-armed patterns present in this simulation. 
The two-armed and three-armed patterns have different ra- 
dial extents and different pattern speeds. Figure [Tip shows 
that the transition between a two-armed and three-armed 
pattern can cause a discontinuity in the spiral arm peaks. 
This type of discontinuity can also give neighborhoods with 
two sets of orbital orientation angles and so gaps in a veloc- 
ity distribution, as illustrated by the blue neighborhood in 
this figure. 

We now compare the list of radial discontinuities and 
gaps with estimated resonance locations identified from the 
pattern speeds in the simulations. If we only had one strong 
wave then we would be restricted to a few resonances. How- 
ever our simulation shows multiple patterns, each with a 
different radial range. Keeping this in mind, perhaps it is 
not surprising that the density distribution shows multiple 
spiral arm discontinuities and the velocity distributions show 
multiple gaps. 

If the spiral patterns are excited at resonances then we 
expect that they will start or/and end near resonances. In 
section 4 we discussed likely ranges for these patterns with 
likely resonances as boundaries. The innermost discontinu- 
ity, seen both as gaps in the velocity distribution and in the 
spiral morphology, lies at radius of about r g ~ 4 kpc. This 
lies near the bar's m = 4 ILR and is not distant from the 
m — 2 ILR of the inner two-armed spiral. The discontinuity 
at r g ~ 5 kpc lies near the bar's corotation resonance, the 
the m = 3 ILR of the three-armed spiral and is not distant 
from the m = 4 ILR of the inner two-armed spiral and the 
m — 2 ILR of the outer two-armed spiral. The discontinuity 
at a radius ra ~ 11 is near the m — 3 OLR of the three- 
armed spiral and the m — 2 OLR of the inner two-armed 
spiral. That at r^S kpc is near the bar's m = 2 OLR and 
the m = 2 ILR of the outer two-armed spiral. It also may be 
the outer boundary of the inner two-armed spiral and near 
its m = 4 OLR. The division at ra, ~ 6.3 is not obviously 
associated with a resonance but could be due to a division 
between the inner and outer two-armed spirals. 

The first order approximation used here has a small di- 
visor problem nearing resonance. To understand the dynam- 



ics near resonance a higher order model is required ( Con- 
topoulos|1975 1. Approaching resonance the epicyclic ampli 



tude increases but does not go to infinity and closed orbits 
exist on both sides of resonance. Consequently a closed orbit 



approximation (such as illustrated in Figure 11 1 may not be 



invalid even quite close to resonance. As epicyclic amplitudes 
become high near resonance, Lindblad resonances are par- 
ticularly likely locations where orbits can cross into nearby 
neighborhoods. Because the orientation of closed orbits can 
vary on either side of a Lindblad resonance they can also 
induce a gap in a local velocity distribution ( Dehnen|2000 |, 
also see Figure 4 by Quillen & Minchev ( 2005 1 . If there are 



two perturbations, localized regions near resonance can be- 



come chaotic ( Quillen 2003 1 perhaps leading to the trun- 



cation of spiral arms. The Hamiltonian model by |Quillen| 
( 2003 1 showed that periodic orbits can exist near resonance 



even when two patterns are present though if there is strong 
resonance overlap there may be regions lacking nearly circu- 
lar orbits where localized diffu sion can be rapid (|Shevchenko| 
2011) and heating is possible (Minchev & Quillen 2006[). 

The gaps in the velocity distribution and discontinuities 
in spiral morphology seen in our simulation can be explained 



with a model containing multiple spiral density waves. The 
onset radii of each wave introduces the discontinuities in the 
spiral arm morphology and associated gaps in the velocity 
distribution. Our simulation shows multiple patterns that 
are coupled, and there are radii that are likely resonant with 
more than one pattern. These radii are likely associated with 
onset radii of waves, gaps in the velocity distribution and 
discontinuities in the spiral arm morphology. 

7.1 Clumps like the Hercules Stream 

When the neighborhood lies on a discontinuity, two sets of 
orbits can pass through the neighborhood; see for example 
the blue neighborhood in Figure [TT| i. In this case the first set 
is those associated the spiral arm peak and their velocities 
would lie on an arc passing through the top right quadrant 
in the uv plane. The other set of orbits comes from the 
inner galaxy and differs in orientation by about 90°. This 
difference in orientation implies that its u values will have 
the opposite sign, and since the orbits have lower guiding 
radii, they should have lower v values. Thus we expect the 
velocities of stars on the second group of orbits can have 
low v and the opposite sign of u and so form a clump on the 
lower left quadrant in the uv plane, similar in position to the 
Hercules stream at v ~ —50 km/s in our galaxy (Dolmen fc 
Binneyp998 l. 

The m = 2 outer Lindblad resonance with the bar at 
earlier times in the simulation is at about 9 kpc. If the solar 
neighborhood lies just outside this resonance with the bar in 
the Galaxy then the radius in this simulation most matching 
our position in the Galaxy (just outside the bar's OLR) is 
at ro ~ 10 kpc. We show velocity distributions in neighbor- 
hoods with this radius and at different times in Figure [l3"| 
At this radius, clumps similar in position to the Hercules 
stream are seen at angles (with respect to the bar) of (j>o = 
90, 135, 270 and 315° but not at all times. It is rare to see 
these streams at angles nearly aligned or perpendicular to 
the bar. This suggests that non-circular motions caused by 
the bar may be needed to produce the large epicyclic ampli- 
tudes corresponding to the strong negative v values in the 
clump. 

The orientation of a periodic orbit can depend upon 
its location. For example, orbits can be aligned with a bar 
outside the m — 2 Outer Lindblad resonance and aligned 
perpendicular to it between this resonance and the corota- 
tion resonance. A flips in orbit orientation on either side of 
a resonance can also cause a gap in the velocity distribution 
(e.g., |Quillen fc Minchev|2005 1. However if the orbital distri- 
bution was primarily determined by the bar then we might 
expect the outer galaxy to look like a ring galaxy. However 
our simulated galaxy at ro ~ 10 kpc has a strong two-armed 
spiral structure. The bar's m = 2 outer Lindblad resonance 
itself may also be a site for coupling between waves. Both 
the bar and local spiral structure affect the Hercules stream 
like clumps in the orbital distribution. 

The velocity distribution in the solar neighborhood not 
only displays a stream at v ~ —50 km/s (the Hercules 
stream) but is triangular shaped for velocities with low u, v 
We show in Figure [14] rightmost panel the velocity distri- 
bution of the magnitude-complete, kinematically unbiased 
sample of 16682 nearby F and G dwarfs by |Holmberg et| 
al. (20091. A triangle shaped distribution at lower u, v and 
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a Hercules stream like feature at negative v are sometimes 
seen in our velocity distributions at ro = 10 kpc, near the 
bar's m — 2 outer Lindblad resonance, are both seen for ex- 
ample the rightmost panel of Figure |10| We have looked 
through our velocity distributions at this ro and with a 
bar orientation consistent with the ~ 45° angle (between 
bar semi-major axis and Sun the Galactic center line) esti- 
mated fr om the SPITZER/GLIMPSE survey ( |Benjamin et| 
al. 2005| l for other examples of velocity distributions that ex- 



hibit the gross features present in the solar neighborhood ve- 
locity distribution. These are shown in Figure [T4| for smaller 
neighborhoods with size set by / = 0.05 along with disk sur- 
face density at the same timesteps in Figure |15| In Figure 
|15| blue circles show the location of neighborhoods used to 
extract the velocity distributions. The disk has been flipped 
so that rotation is clockwise rather than counter clockwise 
and rotated so the extraction neighborhood is oriented at 
the top of the plot. This has been done so that these images 
can more readily be compared with models for the Galaxy. 
The location of the mock Sun is at the blue circle with a 
galactic longitude of 90° viewing to the right. 

At certain times and locations in our simulation the lo- 
cal velocity distribution is similar to that seen in the solar 
neighborhood at a position consistent with the orientation, 
length and distance of the Galactic bar. By similar we mean 
with a triangular shaped feature at low uv and a Hercules- 
like stream at low v. Weak clumps at low u,v are seen in 
our velocity distributions, similar to those labelled as low 
velocity streams such as the Hyades and Pleiades stream 
in the solar neighborhood ( Dehnen fe Binney|1998 1. We do 
not match each of the low velocity streams seen in the solar 
neighborhood velocity distribution but our velocity distri- 
butions do exhibit low velocity clumps suggesting that they 
could be related to spiral waves in the disk. However none 
the models have have discussed above would easily account 
for this structure. Our Hercules like stream is more heav- 
ily populated by stars than that in the solar neighborhood, 
probably because our disk is hotter. Our simulated Hercules 
stream is at lower v than that in the solar neighborhood sug- 
gesting that a somewhat smaller galactocentric radius would 
provide a better analog, however at smaller radii it is more 
difficult to see a strong Hercules like stream in the simulated 
velocity distributions. 

It is interesting to compare the spiral morphologies seen 
in Figure [15] In all cases the location of the mock Sun is 
just exterior to a strong spiral arm. This spiral arm would 
have tangent (as viewed from the mock Sun) at a Galactic 
longitude of about —50° (to the left on the plot), and ap- 
proximately consistent with the location of the Centaurus 



Arm tangent measured by Churchwell et al. ( 2009 1 from the 



GLIMPSE survey counts and denoted as the Crux-Scutum 
Arm tangent by | Vallee~| ([2008 ) . Our simulated galaxy would 
exhibit only one or two peaks or tangent arms in a projected 



number count diagram (similar to figure 14 by Churchwell 
|et al.|[2"009| ) and so would resemble the observed infrared 
observations. As our simulation lacks gas we cannot see ad- 
ditional gas rich but star poor spiral arms that are likely 
present in our Galaxy (Drimmel & Spergcl 2001). While all 



the images in Figure |15| have a strong spiral arm just inte- 
rior to the mock Sun, the opposite side of the galaxy has a 
variety of morphologies. This implies that the solar neigh- 
borhood velocity distribution alone is not sufficient to break 



the degeneracies of a model that includes multiple spiral 
waves. 



8 SUMMARY AND DISCUSSION 

Using an N-body hybrid simulation we have looked for a 
relationship between structure in local velocity distribu- 
tions and spiral and bar density waves. Patterns for these 
waves are identified through Fourier spectrograms. Previ- 
ous studies (e.g., [Rautiainen fc Salo||1999[ ), have suggested 
that there are couplings between waves. The similarity be- 
tween spectrograms at different stages in our simulations 
(as the bar slows down) suggests that our simulation also ex- 
hibits couplings between waves, with bar and lopsided waves 
likely driving spiral waves. Even though multiple waves are 
present, they could be coupled and slowly evolving rather 
than transient. 

We find that arcs in the local velocity distribution are 
seen in neighborhoods located on arm peaks, whereas in- 
terarm regions tend to have lower velocity dispersions. This 
behavior is understood by considering a series of nested ellip- 
tical orbits, appropriate for a cold disk influenced by a single 
spiral perturbation. Spiral arm peaks are present where el- 
lipses touch or overlap and these locations correspond to 
neighborhoods containing orbits with a wider range of ve- 
locity vector orientations than interarm regions. 

Gaps can be seen in the velocity distributions in our 
simulation all over the simulated galaxy, confirming the early 



finding by Fux (20001. However, the v values of these gaps 



correspond to a small number of guiding radii that are near 
resonances with the bar and spiral arm patterns or where 
there are kinks or discontinuities in the spiral arm morphol- 
ogy. The relation between gaps in the velocity distribution 
and discontinuities in spiral structure is understood by con- 
sidering a series of nested ellipses that have a discontinu- 
ity in their relative orientation angles. The discontinuity in- 
troduces neighborhoods containing overlapping orbits that 
have a bifurcated distribution of orbit orientations. This bi- 
furcated distribution corresponds to two clumps and a gap 
separating them in the velocity distribution. Multiple waves, 
each with a different pattern speed, present in the simula- 
tion, introduce radii at which there are discontinuities in the 
spiral structure. Consequently gaps in the velocity distribu- 
tions can be seen at multiple radii in the disk. 

N-body simulations have the advantage over test par- 
ticle simulations that the orbits of particles are consistent 
with the gravitational potential generated by the particles 
themselves. Additionally N-body simulations naturally con- 
tain varying or evolving bar and spiral structures that would 
require additional and possibly poorly constrained parame- 
ters to describe inside test particle simulations. The study 
here, though preliminary and based on only one hybrid sim- 
ulation, can be considered to be complimentary to studies 
carried out with test particle simulations. 

8.1 Implications for Studies of the Milky Way 

From our simulations we have found that interarm regions 
tend to have lower velocity dispersions than arm peaks. They 
also tend to lack arcs and gaps in the velocity distribution. 
The Sun is likely in an interarm region ( Vallee ||2008 1. We 
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may have biased view of the Milky Way disk dynamics be- 
cause the local velocity distribution does not represent that 
at other positions in the Galaxy. For example the mean local 
velocity dispersion in the Galactic disk may have been un- 
derestimated because of our interarm location in the Galaxy. 
Large arcs in the local stellar velocity distribution may be 
present in nearby spiral arm peaks. 

There remains little consensus among different studies 



for the pattern speed of local spiral structure (e.g. Shaviv 
[2003} |Lepine et al.|2010| |Gerhard|2010| . The lack of consen- 
sus may arise because the presence of more one wave com- 
plicates analysis (e.g., Naoz fc Shaviv|2007 \ . The Local Spur 
which is thought to lie between the dominant two-armed or 
four-armed structures may arise because multiple waves are 
passing in and out of phase with each other (e.g., Henry 



et al.|[2003 Meidt et al. 20091. As armlets can appear and 
disappear, a localized short burst of star formation may oc- 
cur that does not continue to progress around the galaxy as 
is expected when a sole wave passes through the disk. Evi- 
dence for short isolated bursts of star formation in the disk 
would support interpretation of the Galactic disk in terms of 
multiple waves simultaneously present in the disk. We have 
noted that density peaks due to constructive interference 
of patterns tend to begin at smaller radius and progress to 
larger radius before fading away. Trails of recently formed 
clusters with age decreasing with increasing radius may be 
found in the Galaxy disk. 

Recent studies have placed strong constraints on the 
pattern speed of the Galactic bar from the morphology of the 



Hercules stream in the local velocity distribution (Dehnen 
|2000} |Minchev et~aT1[2007l |Gardner et"aT1[20Tol ). The as- 
sociation of the v value of the gap separating the Hercules 
stream from the rest of the stars with the bar's m = 2 
outer Lindblad Resonance tightly constrains the bar pat- 
tern speed. However features similar to the Hercules stream 
are seen in our velocity distributions that both associated 
with the bar's m = 2 outer Lindblad resonance and at other 
locations in the galaxy. These clumps do not remain fixed in 
time and are affected by the coupled spiral structure in the 
outer Galaxy. Previous models of the Hercules stream have 
neglected spiral structure and have been based on test parti- 
cle simulations, not N-body simulations (with the exception 
of the study by Fux|2000 1. If wave coupling is common then 
the bar's outer Lindblad resonance may be coincident with 
a Lindblad resonance with a spiral pattern. Future models 
could consider both spiral and bar patterns when interpret- 
ing the velocity distribution near the bar's outer Lindblad 
resonance. 

Velocity distributions with morphology similar to that 
in the solar neighborhood are seen in our simulation at lo- 
cation just exterior to the bar's m = 2 outer Lindblad reso- 
nance with bar orientation consistent with infrared studies. 
By similar we mean, showing a Hercules like stream and 
a triangular shape at lower u, v. The velocity distribution 
more closely resembles that seen in the solar neighborhood 
when there is a strong spiral arm, consistent with the ob- 
served Centaurus Arm tangent, just interior to the solar 
neighborhood. 

Like many other galaxies, the Milky Way is lopsided 
( Levine et al. 2006 ) . If coupling between a lopsided motion 
and a Galactic bar is common then we might expect a strong 
three armed spiral wave would be present in the solar neigh- 



borhood. However, almost all studies have assumed two or 
four armed spiral structures ( Vallee |2008 1. If a strong three 
armed wave is present in the solar neighborhood and its in- 
ner Lindblad resonance is coincident with the Galactic bar's 
corotation radius (as likely true in our simulation) then we 
can estimate its pattern speed. At a inner Lindblad reso- 
nance the pattern speed is related to the angular rotation 
rate by f2 p = f2(l — -%)• If we se t m this relation equal 
to the bar pattern speed (approximate 1.9 times the angular 
rotation rate at the solar circle Qp; |Gardner et al.|2010 l, use 
a flat rotation curve with § = y2 and m = 3 we estimate a 
pattern speed of Q3 = (1 — \/2/3)1.9f20 ~ f2 Q putting the 
Sun near corotation of this three armed pattern. Pattern 
speeds for spiral structure fall into two groups, a faster set, 
placing the Sun near the corotation resonance and a slower 
set placing the Sun near the inner 4:1 Lindblad resonance 



(see Table 3 by Shaviv 20031. A three armed pattern with 
pattern speed similar to f2© would have m = 3 inner Lind- 
blad resonance near the bar's corotation resonance and so 
could be coupled to the bar. In our simulation a two spi- 
ral pattern is also prominent near the position of the mock 
solar neighborhood with inner Lindblad resonance near the 
bar's corotation rotation resonance. We can estimating the 
pattern speed of a two armed wave satisfying this condition 
from the bar's pattern speed using m = 2 finding a pattern 
speed of 0,2 ~ O.55f2o or about 14 km s _1 kpc -1 consistent 
with a slower group of pattern speeds estimated a number 
of s tudies (see Table 3 by |Shaviv|2003| 



Lepine et al. ( 2010 1 propose that a nearby kink in a spi- 



ral arm winding angle is associated with a single Lindblad 
resonance with a spiral pattern. However our simulations 
suggest that kinks and discontinuities in the spiral struc- 
ture not only are associated with Lindblad resonances with 



a single pattern (as also explored by Quillen & Minchev 



2005 



Contopoulos 1988p but are also signatures of transi- 



tions from one pattern to another. Rather than rare, our 
simulation suggests that these discontinuities may exist all 
over the Galaxy. 

The mode or wave coupling possibly present in our sim- 
ulations implies that resonances between different patterns 
could be commonly coincident. As shown by Minc hev fc| 
Quillen ( 2006 1 this causes chaotic motions in the stellar or- 



bital motions and so accelerated heating (or increasing ve- 
locity dispersion) near localized regions of resonance over- 
lap. If waves are coupled in the Milky Way disk then we 
might except both localized accelerated heating and radial 



migration ( |Minchev fc Famaey||2010| |Minchev et al.| 2011 l 
as a consequence (also see |Quillen||2003||Voglis et al.||2006| 



Shevchenko[20TT | 

The waves present in our simulation are likely coupled 
and so not transient. The divisions seen in the velocity dis- 
tribution we interpret in terms of discontinuities associated 
with resonances and transitions from one pattern to another 
rather than a model involving stochastic growth of transient 
spiral structure (e.g., De Simone et al.|2004 i. 



8.2 Future Work 

This study only considered one single simulation and it ran 
only for 1.3 Gyr. The simulation is not ideal as we have un- 
dersampled the bulge and halo in order to well resolve the 
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disk, and the energy conservation is not excellent. However, 
this particular simulation is particularly rich in spiral den- 
sity waves and so a good one for comparing velocity distri- 
butions to structure seen in spectrograms. While our galaxy 
was simulated in isolation, real galaxies suffer external per- 
turbations. Lopsidedness in our simulation that may be an 
artifact of our undersampled live halo, is commonly seen in 
external galaxies (Jog & Combes 2009). Hence the coupling 
between bar, lopsided motion and three armed waves that 
we suspect is present in our simulation may be ubiquitous 
rather than just a curiosity. 



In our simulations we purposely undersampled the halo 
and bulge particle distributions so that we could better re- 
solve the galaxy disk. However as a result our few and rel- 
atively massive halo particles cause spurious heating in the 
disk. We have measured the radial velocity dispersion in the 
disk at a radius of 10 kpc and found that it does not in- 
crease by more than 10 km/s across the entire simulation. 
This suggests that heating by massive halo particles is not 
excessive in this simulation. The similarity between spectro- 
grams in the middle and later part of the simulation suggests 
that waves are coupled to the bar rather than randomly gen- 
erated by noise associated with the granularity of our halo 
during the simulation. Were we to raise the number of mas- 
sive particles in the disk, halo and bulge, the simulations 
would exhibit reduced numerical noise and so might have 
weaker or delayed spiral and bar structures. Future stud- 
ies should not only consider more accurate simulations with 
larger numbers of particles but also explore simulations with 
differently distributed disks, bulges and halos, the role of a 
cold gaseous component and perturbations due to mergers. 



Here we have only qualitatively attempted to relate spi- 
ral morphology and patterns to structure seen in the velocity 
distributions. However future work could aim to bridge the 
gap between test particle simulations and N-body simula- 
tions by understanding the orbits in each feature seen in the 
velocity distributions. Here we have focused on gross fea- 
tures present seen in the velocity distributions, such as gaps 
and arcs. Future studies can also aim to understand the finer 
features as well, that are primarily seen in smaller neighbor- 
hoods (e.g., see Figure 10 1. These studies would help in- 



terpret forthcoming pencil beam (e.g., Minchev & Quillen 
2008 ) and large scale velocity and proper motion surveys of 



stars in the Galactic disk. 



Here we have interpreted gaps in the velocity distri- 
butions in terms of changes in the shapes of nearly closed 
orbits associated with different density waves. However fu- 
ture studies could also consider alternate interpretations for 
structure seen in velocity distributions such as transient spi- 
ral arms ( |De Simone et a l. 2004), bifurcations in the fami- 
lies of periodic orbits ( Katsanikas et al.|2011 1, chaotic orbits 
along spiral arms near a bar ( Voglis et al.|2006 1 , phase wrap- 
ping following an initial disturba nce (jMinchev et al.|2009[) or 
individual Lindblad resonances flQuiilen fe Minchev|20"05l ). 
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Figure 1. Disk stellar densities projected into the xy plane, during early spiral arm and bar growth. The time of each snapshot is shown 
in Myr on the lower right of each panel, a) Here we show the earlier part of the simulation. During this time we see initial stochastic 
spiral arm growth followed by the onset of bar instability, b) Here we show snapshots from the later part of the simulation. We can see 
that the bar is not symmetrical; it can be lopsided. At times the galaxy appears to have three-arms (e.g., at T = 750 Myr). Kinks and 
discontinuities are commonly seen in the spiral structure, c) The same as b) except only massive disk particles are shown. There is no 
significant difference between the distributions of massive and massless disk particles. 
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Figure 2. Differential surface density as a function of the log of the radius r and azimuthal angle <j>. Here we show the surface density 
at three different times early in the simulation and during spiral arm growth. The times of the snapshots in Myr are shown in the upper 
right hand corners of each plot. Patterns move to the right with time and trailing structures have negative slopes. A logarithmic spiral 
arm has high densities along a line on this plot with pitch angle given by the slope of the line. The vertical feature in the leftmost panel 
tilts and moves to the right in the center and then rightmost panel. In the rightmost panel it is a stronger perturbation and with a larger 
and trailing pitch angle. The spiral arm shown in the center of the leftmost panel becomes more tightly wound, and stronger, as would 
be expected from a swing amplification mechanism. 
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Figure 3. Similar to Figure [2] except showing the differential surface density after bar formation. On the lower parts of these panels the 
vertical features are the bar. There are two vertical lines per panel as would be expected from an m = 2 bar-like structure. Exterior to 
the bar, both two-armed and three-armed patterns are seen and their patterns move more slowly than the bar. On an animated version 
of this plot it is clear that patterns in the inner galaxy are faster than those in the outer galaxy. Interference between the patterns 
sometimes causes armlets to appear and disappear. 
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Figure 4. Spectrogram of the m = 2 Fourier components in a window spanning a) from time T = 350 to 950 Myr and b) from time 
T = 750 to 1305 Myr. The y axis is angular frequency in Myr -1 and the x axis is radius in kpc but shown on a log 10 scale. Overplotted 
are five lines that are twice the following angular frequencies in order of bottom to top fi — k/2, Q — k/4, f2, Q + k/4, and Q + k/2. 
The bar, once present, is stable and has a slowly decreasing pattern speed. There are at least two two-armed spiral patterns. The spiral 
patterns at both early and later times seem to extend between their inner m = 2 or 4 Lindblad resonance out to their corotation radii. 
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Figure 5. Spectrogram of the m = 4 Fourier component in a window spanning from a) time T = 350 to 950 Myr and b) from 750 to 
1305 Myr. Similar to Figure|4] With the exception of the bar (with u = 0.016 at early times and 0.013 Myr -1 at later times) the spiral 
patterns seen here are strongest between their inner and outer m = 4 Lindblad resonances (as seen via intersections with green lines 
showing Aui ± k. Higher frequency features are seen above the bar's angular frequency at both early and later times. 
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Figure 6. Spectrogram of the m = 3 Fourier components in a window spanning from a) time T = 350 to 950 Myr and b) from 750 to 
1305 Myr. Similar to Figure [1] but constructed using the m = 3 Fourier components instead of the m = 2 ones. Overplotted are seven 
lines that are three times the following angular frequencies in order of bottom to top f! — re/2, f! — re/3, Q — re/4, Q, Q + re/4, f! + re/3, and 
f2 + re/2. The dominant three-armed pattern extends between its m = 3 inner Lindblad resonance and its m = 3 outer Lindblad resonance 
at both early and later times. The bar itself on this plot has an angular frequency of u)^ f» 0.120 at early times and uif, f» 0.105 Myr -1 
at later times. There is a possible coupling between the corotation resonance of the bar and the m = 3 inner Lindblad resonance of the 
three-armed spiral pattern as these two resonances are close. 
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Figure 7. Spectrogram of the m = 1 Fourier component in a window spanning from a) time T = 350 to 950 Myr and b) from 750 to 
1305 Myr. Similar to Figure|6] The bulge centroid position was not subtracted previous to measurement of this spectrogram. On the left 
the bright feature at u) = 0.08 Myr -1 is associated with the bar. The bar can be seen in the m = 1 spectrogram because it moves with 
the lopsided oscillation. Low frequency features may be symptomatic of non-linear coupling between waves. 
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Figure 8. Local velocity distributions at different positions in the outer galaxy at a simulation time T = 600 Myr. Velocity distributions 
are extracted from neighborhoods centered at 8 equidistant angles and at 6 radii spaced logarithmically. Each panel shows the distribution 
in u (x-axis) versus v (y-axis). Each row shows neighborhood at a particular galactocentric radius with the lowest row with ro = 4.1 kpc 
and the top row with ro = 12.5 kpc. Each column shows neighborhoods with a particular galactocentric angle with the leftmost aligned 
with the bar. Each neighborhood is separated in angle by 45° with angle increasing counter-clockwise in the galaxy plane. Positions and 
sizes of each neighborhood are shown in Figure [9] Here we see that clumps exist in the velocity distribution at all radii in the galaxy. 
Features in the velocity distribution at one neighborhood are often seen in the nearby neighborhood at a larger or smaller radius but at 
a shifted v velocity. 
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Figure 9. The positions of the neighborhoods used to construct velocity distribution histograms shown in Figure H] a) The neighborhoods 
are shown in a Cartesian coordinate system. The green circles show the angle of the neighborhoods oriented along the bar and referred 
to as angle zero (with respect to the bar) and are shown in the leftmost column of Figure [8] In Figure [8] the angles of each column, 
from left to right, correspond to positions with angle increasing by 45° in the anti-clockwise direction on this plot. Radii of changes 
in spiral arm pitch angle or discontinuities in spiral structure estimated from the polar plots are shown as dotted black circles, b) The 
approximate locations of the neighborhoods shown on a projection of angle versus logarithmic radius (the differential density distribution 
is constructed as in Figure |3j. The positions of the neighborhoods on this plot are in the same order as the velocity distribution panels 
shown in Figure [8] Black dotted horizontal lines are shown at radii where there are changes in spiral arm pitch angle or there are 
discontinuities in the spiral arms. Many of these lines correspond to resonances with patterns seen in the spectrogram. 




Figure 10. A comparison between velocity distributions extracted at a radius of rg = 10 kpc using neighborhoods with radius with 
500 pc (/ = 0.05; top panels) and 1000 pc (/ = 0.1; middle row of panels). The simulation time was arbitrarily chosen and is 800 Myr. 
The bottom row is similar to the middle panel but shows only massive disk particles whereas the distributions in the top two rows are 
constructed from both massive and tracer particles. The neighborhoods have the same angles as in Figure [8] and are arranged in the 
same way. When the neighborhoods are larger the velocity distributions are smoother and less grainy as would be expected because more 
particles are found in each neighborhood. 
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Figure 11. a) A series of ellipses is shown with semi- major axes logarithmically spaced from the other. For a cold disk, to first order in 
the perturbation strength and for an m = 2 perturbation, orbits are approximately ellipses. The orientation of each ellipse semi-major 
axis is 20 degrees counter-clockwise lower than that just inside it. The only exception to this rule is between the 6-th and 7-th ellipse 
(counting from smallest to largest) where the angular difference is larger. Spiral arm peaks (density peaks in a stellar distribution) occur 
where the orbits are close together or overlap. The red circle is an example of an arm peak. There are a large range of velocities or orbit 
orientations in a region where many elliptical orbits overlap, giving an arc in the velocity distribution. In interarm regions, (e.g., green 
circle) the range of velocities (or angles on this plot) in a particular neighborhood is low and so the corresponding velocity dispersion 
in this neighborhood is also low. The spiral arms are trailing with galactic rotation in the counter clockwise direction. On an arm a 
broad range of angles exists in the orbits crossing the neighborhood, so the arc in the corresponding velocity space is continuous. The 
angular offset between the 6-th and 7-th ellipse, as might occur if there are different patterns in different regions of the disk, causes a 
discontinuity in the overlap region that manifests as a gap in the velocity distribution, b) The inner 6 curves are ellipses. The outer 
ones are a sum of elliptical and triangle perturbations such as might occur if there are simultaneously two-armed and three-armed spiral 
density waves in different regions of the disk. The transition region exhibits discontinuities in overlap regions caused by the onset of 
the three-armed pattern. The blue circle is located at a radius where there is a discontinuity and shows that the angular distribution of 
orbits is bifurcated. We expect that in such a region there will be gap in the velocity distribution. 
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Figure 12. We show uv velocity distribution histograms generated with neighborhood radii of ro = 6.4 kpc but at different times. 
From top to bottom row the histograms were generated from snapshots at output times 500 to 950 Myrs. Each snapshot is separated 
in time by 50 Myr. As in Figure [8] the angles are given with respect to the bar and the same sequence of angles is shown. Gaps in the 
velocity distribution tend to remain at the same v values. Even though the positions shown are fixed in the bar's rotating frame there 
are variations in the features seen in the velocity distribution that must be due to the spiral waves in the outer disk. 
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Figure 13. Velocity distribution histograms generated with neighborhood radii of ro = 10.0 kpc but at different times. Similar to 
figure [l2] These neighborhoods lie just outside the bar's m = 2 outer Lindblad resonance (at about 9 kpc). A clump similar to the solar 
neighborhood's Hercules stream (with both negative u and v) is primarily seen at angles <j>o = 90, 135, 270 and 315° respect to the bar. 
It is rare to see these streams at other angles with respect to the bar. Non-circular motions caused by the bar may be required to induce 
the large epicyclic amplitude implied by the low v values in the clump. 
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Figure 14. The leftmost 4 panels show velocity distributions at a radius of tq = 10 kpc chosen because they exhibit both a Hercules-like 
stream and a triangular shape near the origin similar to the morphology of the solar neighborhood's velocity distribution (shown on the 
right). These velocity distributions where created with neighborhoods with / = 0.05. The time (in Myr) of each snapshot is on the lower 
left of each panel. 




-10 -5 5 10 -10 -5 5 10 -10 -5 5 10 -10 -5 5 10 
x (kpc) x (kpc) x (kpc) x (kpc) 



Figure 15. For the four velocity distributions shown in Figure 1 1 4| we show the morphology of the disk and the location (as small 
blue circles) of the neighborhoods used to create these distributions. The disk images have been flipped so that rotation is clockwise 
and rotated so that extraction neighborhoods are on the top. This has been done so that these figures can more easily be compared to 
maps of the Galaxy. The blue circle represents a neighborhood chosen to be like the solar neighborhood. As our simulated bar is longer 
than the Galactic one, the chosen neighborhood is somewhat more distant than the Sun from the Galactic center. The bar orientation is 
consistent with that inferred from infrared observations of our Galaxy. These images have in common a strong spiral image just interior 
to the mock solar neighborhood. This arm would be consistent with the Centaurus Arm tangent arm. 
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Table 1. Estimated Pattern Speeds and Locations of Resonances 



Pattern 


Time 


Pattern speed 


Tin Tout 


ILR 2 


ILR 3 


ILR4 


CR 


OLR4 


OLR3 


OLR 2 




(Myr) 


(radians/Myr) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


Bar 


350-950 


0.040 


0.5 -6 




2.4 


3.6 


5.1 


7.2 


7.6 


9 


Inner two-armed Spiral 


350-950 


0.030 


3 -10 


2.4 


3.8 


1.4 


7.2 


9.5 


10 


11 


Outer two-armed Spiral 


350-950 


0.015 


6 -15 


4.7 


7.1 


9.1 


13.5 


17 


18 


19 


Three-armed Spiral 


350-950 


0.023 


3 -15 




1.5 




9.1 




12.3 




Bar 


750-1305 


0.035 


0.5 -6 




3.6 


3.0 


6.0 


8.1 


9.1 


10 


Inner two-armed spiral 


750-1305 


0.022 


3 -10 


3.2 


1.9 


6.0 


9.5 


12.0 


13.3 


14 


Outer two-armed Spiral 


750-1305 


0.013 


6 -20 


6.3 


8.9 


10.9 


15 


19 


19.5 


22 


Three armed Spiral 


750-1305 


0.023 


5 -15 




4.5 




9.1 




12.3 





These pattern speeds are measured from features seen in the spectrograms shown in Figures [4] and [6] The pattern speeds can be 
approximately converted to units of km s _1 kpc -1 by multiplying by 1000. The times listed in the second column correspond to the 
range of times used to construct these spectrograms. The radii r^ n and r ou t are estimated inner and outer radii in kpc of the individual 
patterns. Here ILR2 and OLR2 refer to the inner and outer m = 2 Lindblad resonances. Likewise ILR3 and ILR4 refer to the the m = 3 
and m = 4 inner Lindblad resonances. When a resonance from one pattern is near that of another pattern, one pattern could be 

coupled to or driving the other. 



